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Abstract. A second-order differential (q- difference) eigenvalue equation is 
constructed whose solutions are generating functions of the dual (q-)Hahn 
polynomials. The fact is noticed that these generating functions are reduced 
to the (little q-) Jacobi polynomials, and implications of this for quasi-exactly 
solvable problems are studied. A connection with the Azbel-Hofstadter prob- 
lem is indicated. 
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1 Introduction. 



In the present paper we shall consider some second-order differential and 
g-difference eigenvalue equations of the form: 

a(z)f"(z) + b(z)f'(z) + c{z)f{z) = \f(z) (1) 

and 

a{z)f{q s z) + P(z)f(q s+1 z) + ^z)f(q s+2 z) = \f{z) (2) 

respectively, where the functions a(z), b(z), c(z) are polynomials in z, while 
a(z), /3(z), j(z), in z and z~ l . We shall be looking for polynomial solutions 
f{z). Note that using a transformation of the type ip(y) = g(y)f{z(y)) we 
can always reduce (1) to the Schrodinger form — -^ip{y) + V{y)il){y) = \ip{y)- 

After subjecting the coefficients of (1) and (2) to certain general condi- 
tions, we shall see that the polynomial solutions to equations (1) and (2) are 
given by generating functions of the dual Hahn and dual g-Hahn polynomials, 
respectively. These solutions are explicit in the sense that all the eigenval- 
ues A and corresponding eigenfunctions f(z) in the space of polynomials of 
degree at most iV are known explicitly. 

Let us first consider equation (1). It is known |IJ (see also 0] for a 
recent review of the subject) that the spectral problem for the operator D = 
a ^ z )'^ Jr ^ z )'h JrC ^ z ) m ^ ne space Hn spanned by the vectors 1, z, z 2 , . . . , z N 
is closely related to the representation theory of the algebra sl 2 - 7~Cn is a 
representation space of this algebra, and the generators of sl 2 have in this 
representation the following form: 

J + = z 2 ^-Nz- J° = z^--N/2; J~ = A ( 3 ) 

dz dz dz 

The necessary and sufficient condition || for the operator D to leave TCn 
invariant is that D be expressed in the form 

D = c ++ J+ J+ + c +0 J + J° + c+_ J + ,r + c _ J°J- + 
c—J~J~ + c + J + + c J° + c_J~ +d, 



where Cij,Ci,d are constant parameters. Henceforth, we shall assume that 
(4) holds. According to the classification of Turbiner, equation (1) is called 
in this case quasi-exactly-solvable.[] 

1 For the expressions of all possible potentials in the corresponding Schrodinger equa- 
tions see j|. 
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If the parameters c +H 



c+o 



c + = then .D is, obviously, upper 



diagonal in the basis of monomials {z k }^ =0 , and hence, it preserves the flag 
TLq C Tti C • ■ ■ C Hn- The coefficient a(z) in this case is a polynomial of no 
more than the second degree, b(z), the first degree, and c is independent of 
z. Hence, the operator D preserves for any N. Corresponding equation 
(1) is then called exactly-solvable^. It is easy to verify that in this case, 

changing the six remaining parameters c + _, c _, c , c , c_, d, we can obtain 

for any iV an arbitrary operator D with the just mentioned restriction on 
the degrees of a(z), b(z), and c(z). The full classification of the polynomial 
solutions to (1) for such an operator D is available in the literature (e.g., 
U). All these solutions can be written in an explicit form. In particular, 
the classical orthogonal polynomials (Jacobi, Laguerre, and Hermite) satisfy 
exactly solvable equations of the type (1). In the present paper, however, we 
shall be interested in a different type of solutions. 

Take once again the general case (4) and let now c ++ = c = 0. Then 

it is seen from (3) and (4) that the operator D is represented in the basis of 
monomials {z k } k=0 by a tridiagonal matrix. The equation (1) takes on the 
following form (f(z) = EfcU^VOQ 



/ a - X b 

Ci ai — X b\ 

c 2 a 2 — A b 2 



\ fPo\ 
Pi 







(5) 



V c N a N -XJ \p N J 

We see that the quantities pk satisfy the three-term recurrence relation 
c fc p fc _i + (a k - X)p k + b k p k+1 = 0, p_i = 0. 



(6) 



Thus [e.g., |?J, Pk(X) form a finite system of orthogonal polynomials. (See 
[0)BBE9 f° r studies related to this aspect of quasi-exact solvability.) 

The matrix elements a k , b k , c k are polynomials of degree 2 in the index 
k. The coefficients of these polynomials are expressed in terms of the 7 free 
parameters c + o, c + _, cq_, c+, Cq, c_, d. It is easy to verify that a k , b k , c k can 



2 An important algebraic approach to the exactly solvable problems is proposed in 0] . 

3 The quantity J^fc=o a k zk Pk, where aj~ are some parameters, is called a generating 
function of the sequence {Pk} k =v 
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be obtained in the following way. Assume A(k), B(k), C(k) to be arbitrary 
polynomials in k of degree 2 and impose the boundary conditions C(N+1) = 
B(-l) = 0. Then a k = A(k), b k = B(k), c k = C(k), k = 0,...,N. 

Now consider a particular case of (5). Namely, impose the Askey- Wilson 
condition for the transposed matrix a k + b k -i + c k+ \ = and the restriction 
c +0 = c _. Then the polynomials {j3 fc (A/c + o)}£L associated with the trans- 
posed matrix, are the dual Hahn hypergeometric orthogonal polynomials. 
One of their generating functions provides an explicit solution to (1). In Sec- 
tion 2 we show that the corresponding equation (1) is reduced to the exactly 
solvable (see above) equation for the Jacobi polynomials by means of the 
transformation ip(y) = (y + 1) N f (^^■) . (Note that this is a particular case 
of the transformation that connects various forms of a quasi-exactly solvable 
equation [|1(J. Such transformations comprise the irreducible representation 
of the group GL{2) in the space of polynomials.) Equation (1) in this case 
has an infinite number of formal explicit solutions f(z), but only N + 1 of 
them are guaranteed to be polynomials. 

The above considerations can be generalized for equation (2). Equation 
(2) is related to quantum deformations of the sl% algebra in a similar way as 
equation (1) is related to sh [p,|TB| Q. For connections between (quantum) 



groups and orthogonal polynomials see 0,0]- 

Let a(z), f3(z), j(z) in (2) be the first order polynomials in z and z~ x . 
Then, obviously, in the basis of monomials {z k } k y =0 , equation (2) has the 
form (5) where a k , b k , c k are expressed in terms of the coefficients of the 
polynomials a(z), (3(z), j(z). For the space spanned by {z k } k=Q to be an 
eigenspace, these coefficients should be such that 6_i = c^+i = 0. This 
condition implies that only 7 out of total 9 coefficients are independent. We 
shall fix 3 more by requiring the Askey- Wilson condition for the transposed 
matrix to hold: a k + b k -\ + c k+ \ = 0. Fixing then one of the remaining 4 free 
parameters in an appropriate way and putting s = 0, we obtain the dual q- 
Hahn basic hypergeometric polynomials {p k (\/e)} k=0 as a system generated 
by the recurrence relation + (a k — X/e)p k + c k+ if) k+ i = 0. 

In Section 3 we shall consider equations of the type (2) whose solutions 
will be given by the generating functions of the dual g-Hahn polynomials. 



4 In addition to the results for differential and g-difference equations reviewed in B], 
it is also possible [0 to obtain similar results for difference equations of the form 
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These equations are related to the g-difference equation for the little g-Jacobi 
polynomials. 

It is interesting to note that the zeros of polynomial solutions of equations 
(1) and (2) are connected with the eigenvalues A by a set of Bethe-ansatz 



type algebraic equations [^5|16 



Thus, the message of the present paper can be summarized as follows. 
The dual Hahn (dual g-Hahn) polynomials are the most general system 
in the Askey-scheme of the known hypergeometric (basic hypergeometric) 
orthogonal polynomials [|14] whose generating function provides an explicit 
polynomial solution to the eigenvalue equation (1) (equation (2)). This gen- 
erating function of the dual Hahn (dual g-Hahn) polynomials is reduced to 
the Jacobi (little g-Jacobi) polynomials. (The corresponding problems are, 
thus, exactly solvable.) 



2 Dual Hahn polynomials and a differential 
equation 

Unlike in the introduction, we shall now begin with the dual Hahn polyno- 
mials rather than with the differential equation. 

The dual Hahn polynomials are defined by the formula (e.g., [II]) 

where 7 and S are fixed parameters and the "shifted" factorial is defined as 
(a) = 1, (a)fc = a(a + 1) • • • (a + k — 1), k = 1, 2, . . .. The polynomials 
(|7|) satisfy the recurrence relation (which we will formally consider for an 
arbitrary integer n) 

X(x)p n = A n p n+1 

~ (A n + C n )p n + C n p n _i, 
A n = ( n - jV)(n + 7 + l), C n = n(n - 5 - N - 1), (8) 
A(x) = x(x + 7 + 5+ 1). 

The above three-term recurrence relation can be viewed as the eigen- 
value equation for an infinite tridiagonal matrix M, p„'s being components 
of an eigenvector. For what follows, we would need to demand that the fi- 
nite dimensional space Cn corresponding to the indices n = 0,1, ... ,N be 
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invariant under the action of the matrix M. This would be the case if the 
matrix elements M_io = Mn+in — 0. Since for our matrix Mo-i = Co = 
and M NN+ i = A N = 0, the transposed matrix M T will have the desired 
property of preserving Cn- The polynomials associated with M T satisfy the 
recurrence 



K x )Pn = C n+1 p n+1 - (A n + C n )p n + v4 n _!p n _i; (9) 
and it is easy to show by induction that 

A Q A l ■ ■ ■ A n _x (-AQ n ( 7 + l) w 

Pn — ^ ^ „ Pn — , r ,,v . Pn- 

CxC 2 ---C n (-5 - N) n n\ 

Now multiply both sides of (§) by z n and perform summation over n from 
n = to N. We obtain 

Xf(z) = z(z-l) 2 f"(z) + 

{( 7 - iV + 2)^ 2 - ( 7 - 5 - 2N + 2)z - (5 + N)}f'(z)- (10) 
iV( 7 + l)(z-l)/(z), 

where f(z) = J2n=o zn Pn- To get the homogenious equation fllPp, it was 
necessary to put p_i = and p^v+i = (we can do this because we are 
looking for solutions in Cn). 

We can represent ([TOj) in the form Xf(z) = Df(z) as the eigenvalue 
equation for a second-order differential operator D in the space Ti^ spanned 
by monomials {z k }^ =0 . One can already notice that fllQ|) can be reduced to a 
hypergeometric equation. However, we shall follow another approach which 
can be easier generalized to g-difference equations. 

Since M T in Cn is just the matrix representation of the operator D in 
the basis {z k }^ =0 , the eigenvalues of D in Hn and M T in Cn are the same. 
To find them, first replace the parameter iV in (|7]) and (|8|) by N + e, e ^ 0. 
Then (^) will be valid not only for n = 0, 1, . . . , N, but also for n = N + 1. 
We find the eigenvalues from the equation: 

= det(M T - XI) = det(M - XI) = \im A Ax ■ ■ ■ A N p N+1 (X). (11) 

Here (only) one of the factors Ai goes to zero as e — ► 0: A N = — e(iV + 7 + 1). 
Furthermore, only the addend with the index k = N + 1 in the expression 

„ (XM) N ^ (-N-l) k (-x) k (x + 1 + 5+l) k 

P N+1 {A[X)) - 7 — TV7 r~Tj 

t^o (7 + l)k(-N - e) k k\ 
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is not bounded as e — > (growing as 1/e). Hence (|Tl]) is equivalent to 
(— x)at+i(x + 7 + 5 + l)jv+i = 0. From here, using the definition of A in (J8|), 
we obtain the eigenvalues: 

A(m) = m(m + 7 + 5 + 1), m = 0, 1, . . . , N. (12) 

The corresponding eigenvectors of D are 

Notice that this generating function is one of those admitting representation 



in terms of the hypergeometric series [14 



/■w-(i-.rs ( "-ff^ +1 V m 

Using one of the representations of the Jacobi polynomials (see, e.g., ||17|| ) 
we can rewrite (|i~3|) in the form: 

(N-m)\ n . N ( _ 5 _jv-i ,-7-jv-i) ( 1 + z 



(-iV-d)jv-m \1-Zj 

where P^ a '^(x) is the usual notation for the Jacobi polynomial. Thus, this 
formula expresses the generating function of the dual Hahn polynomials in 
terms of the Jacobi polynomials. It is easy to transform (|TUp into the hyper- 
geometric equation for the Jacobi polynomials. Since the equation for the Ja- 
cobi polynomials is valid for an arbitrary large index k, formula flU) provides 
an infinite number of nonpolynomial solutions to ( pTOj ) for m — — 1, —2, .... 

The operator D is expressed in the following form in terms of the gener- 
ators (4): 



D = J + J° - 2J + J- + J°J- + (7 + 1 + N/2) J + + (5 - 7 - 2) J°- 
(N/2 + 5)J~ + iV(<5 + 7)/2 

After the transformation 



(15) 



ip(y) = /(coth 2 -) | yj sinh y sinh 7 — cosh 5 - coth^ — 1 (16) 
2 L 2 2 2 J 
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equation ([H]) is reduced to the Schrodinger-type equation 
d 2 



dy 



ip(y)+V(y)ip(y) =eip(y), 



V(y) = .\ 2 {(7 - S)(2N + 7 + S + 2) coshy+ (17) 
2 sinn y 

(JV + 7 ) 2 + (iV + 5) 2 + 2(2iV + 7 + 5) + ^| + l(l + 7 + S) 2 , 
with formal solutions: 



^m(y) 



-m(m + 7 + (5 + l), m = . . . , -1, 0, 1, . . . , JV, 

= ^-^"^^"^(-coshy) (18) 

V / smhy sinh 7 | cosh 5 | coth^ |(1 — coshy) 



A? 



where c m is a constant factor. This is the exactly solvable Schrodinger 
equation related to the Jacobi polynomials. It is easy to verify that if 
7 + JV < and 2m + 7 + 5 + 1 > 0, then the function ipm(y) belongs to 
the space L 2 (— oo,oo) (that is ff^, \ip m \ 2 dy < oo). In this case, since the 
operator —d 2 /dy 2 + V{y) is symmetric with respect to the inner product 
(f-,9) — f^°oc f(y)9(y)dy, such eigenfunctions ip m (y) corresponding to differ- 
ent e m are orthogonal. 

Note that, generally, ip m (y) and ip' m (y) are discontinuous at y = 0. Con- 
sider the physically more reasonable Schrodinger equation (|17D in the space 
L 2 (0, oo) with the boundary condition ip(0) = 0. The functions ip m (y) for 
m = 0, 1, . . . , JV belong to this space if 7 + JV < -1/2, 7 + 5 + 1 > 0. The 
corresponding e m are the levels of the discrete spectrum because they are 
less then the asymptotics (l + 7 + <5) 2 /4 of the potential as y — > 00; and V(y) 
goes to —00 as y — ► 0. Note that e m , m = JV, . . . , 1, are the first JV + 1 
lowest eigenvalues of the Schrodinger operator. 

3 Continuous dual g-Hahn (dual g-Hahn) 
polynomials and a (/-difference equation 
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3.1 q a. root of unity 

The continuous dual g-Hahn polynomials (which depend on parameters a, b, 
and c) are defined by the expression (e.g., [|i4f ) 

Pn{x) = 2^ — rrn — w — x — , 2x = t + t , n = o, 1,..., 

^ {ab;q) k {ac;q) k {q;q) k 

(19) 

where (d; g) = 1 and (d; q) k = (1 — — dg) •••(! — <ig fc A; = 1,2,.... 
(In fact the n'th continuous dual g-Hahn polynomial differs from p n (x) by a 
constant.) They satisfy the recurrence relation 

2xp n = A n p n+1 + (a + cT 1 - A n - C n )p n + C„p n _i, 

A n = a' 1 (l-a6g")(l-acg n ), C n = a(l - g n )(l - bcq" 1 ' 1 ). ~ 

As in the previous section, introduce the matrix M 9 associated with the 
eigenvalue problem (p0|) and the space corresponding to the indices n = 
0,1, ... ,k. Let us first consider the case when q is an iV'th primitive root 
of unityQ, that is q = e 2mS ^ N , where S and N are positive integers which do 
not have a common divisor other than 1. Let us set furthermore ac = q. 
Then, obviously, M q preserves Cn-i- (Moreover, the orthogonal complement 
of £at_i to the whole space where M q acts is also invariant with respect to 
M q .) 

Multiplying both sides of the recurrence relation ( p0|) by z n and perform- 
ing summation from n = to iV — 1, we obtain 

2xf(z) = {(az)- 1 + az}f(z) + 

{-(a' 1 + bq-^z' 1 + a + 2b + qa' 1 - (a + b)qz}f{qz)+ (21) 
b{{qz)- l -q-l + q 2 z}f{q 2 z), 

where /(z)=En=o^>n. 

Proceeding in a similar way as in Section 2, we obtain the following set 
of solutions to (|21 ) in the space spanned by {z k }^~Q-. 

2x m = aq m + a~ l q- m , (22) 



Note that the basic hypergeometric polynomials for q a root of unity have a number 
of interesting properties and applications pq] . 
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f ( \ _ v^ 1 » v (g" n ;g)fc(a 2 g m ;g)fc(g" m ;g)fcg fc 
/ml " j „V fro K ?)*((?; g)*) 2 

(qz; q) N - 1 - m L ^ ^ > m = 0, 1, . . . , JV - 1, 

fc=o {ao,q)k{q,q)k 

(23) 

where in the last formula we used the equivalence of the continuous dual 
g-Hahn polynomials at ac = q~ N+1 and the dual g-Hahn polynomials (to 
be verified below) and the expression for a generating function of the dual 
g-Hahn polynomials jl4] . 



The solution is especially simple for m = 0: 2x = a + a -1 , fo{z) = 
1 + z + z 2 + • ■ • + 2; 7V_1 . In this case we also know explicitly the zeros of fo{z) '■ 
Zi = g\ i = 1, 2, . . . , N — 1. Note that the zeros of all N solutions f m (z) can 
be found in the case when b = 0. Then it is a simple exercise to obtain, using 
(j2~iP , the set of zeros z(m) = {zi, z 2 , ■ ■ ■ , z^-i} of f m (z): 

z(m) = {q m+ \ g m + 2 , . . . , q N ~\ a~ V m+ \ a~ V m+2 , . . . , a" 2 }, 

m = 1,2, 3,..., TV -2 (24) 

z(0) = {g, g 2 , . . . , g^" 1 }, *(JV - 1) = {a^ 2 g 2 , a" 2 g 3 , . . . , a~ 2 q N } 

The difference operator D q (defined by the equation ([H]) written in the 
form 2xf(z) = D q f(z)) can be expressed in terms of the generators of the 
[71/2(5/2) algebra represented in Hn-i- 111 a certain representation in this 
space the generators have the form (we use the notation from [16]): 

N-l JV-1 

A = q~—T + , D = q—T_, 

B = ztfl* - g- 1 / 2 )- 1 (g - <T^T + ), (25) 
C = -^ 1 (g 1/2 -g- 1/2 )- 1 (T_-T + ), 

(recall that g = e 2mS ^ N ) where the operators T + and T_ act on a vector 
g(z) G Hn-i as follows: T±g(z) = g(q ±l l 2 z). 
As is easy to verify, 

A, = A 2 {-b(l + g^ 1 )^ + (g 1 / 2 - g-V2)( 6g - V-iCA + aq'^^BD- 
bq^BA - aT 1 q !i z L+1 CD) + (a + 2b + a' 1 q)q E ^ 1 }. 

(26) 
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3.1.1 Azbel-Hofstadter problem 

It was recently shown |19|] that part of the spectrum of the Hamiltonian in 
the Azbel-Hofstadter problem (of an electron on a square lattice subject to 
a perpendicular uniform magnetic field) can be obtained as (N) solutions A 
of the following equation in TCn-i- 

liz- 1 + qz)f(qz) - liz- 1 + <tM/(<TM = A/(z), (27) 

where q = e i$/2 . $ = 4irS/N is the flux of the magnetic field per plaquette 
of the lattice. (Henceforth, we assume that iV is odd.) The spectrum has 
particularly interesting properties when S, N — > oo so that S/N — > a, where 
a is an irrational number (see, e.g., [2"0|,2~I|). Representation of (WJt) in the 
basis of monomials gives: 

t ( q n+l_ q -(n+l) )pn+i+t ( q n_ q -n )pn _ i = ^ n = , 1 , . . . , iV - 1 , (28) 

where the polynomials p n (X) are defined by the formula f(z) = 

Y,n=0 Z n p n (X). 

On the other hand, setting in ( p0| ) a = iq 1 ^ 2 (hence, c = —iq 1 ^ 2 ), 6 = 0, 
we reduce (|20"1) to 

(1 - q n+1 )p n+ i + (1 - q n )p n ~i = 2xp n , ^ 
p n = a~ n p n , n = 0, 1, . . . , N - 1, 

If we denote the N x N matrices corresponding to eigenvalue equations 
( |28|) and (|29| ) by H and M, respectively, then the following expression holds: 

H = (M — M*)/i. (30) 

In other words, H is the imaginary part of 2M. (Note that M and its adjoint 
M* do not commute.) The spectrum of M is given by ( ^2[) with a = iq 1 ! 2 : 
2x k = 2 sin k = 0, 1, . . . , JV - 1. 

Expression ( ^0|) provides a connection between the results of Section 3 
and the Azbel-Hofstadter problem. 



3.2 Arbitrary q 

Equations (|1^) and (|20|) are valid for an arbitrary complex q (except for 
certain fixed values which one can treat on the basis of continuity consider- 
ations). In this general case, in order to obtain a g-difference equation with 
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the largest number of free parameters, we shall use the approach of Section 
2. Namely, consider the polynomials associated with the transposed matrix 
Mj . Put ac = q 1_N , then the space £jv-i w iU be invariant with respect 
to M g T . (Note that unlike for q N = 1, in the general case the orthogonal 
complement of Cn-i to the whole infinite-dimensional space is not invariant 
with respect to Mj.) The polynomials associated with Mj are connected 
with the dual g-Hahn polynomials as follows (c.f. Section 2): 

(06; q)n(q~ N+1 ;q)n 
Pn a^( q] q) n (ba-iq- N ^;q) n Pn - 

Proceeding as in Section 2, we obtain the following equation for the gener- 
ating function f{z) = T,n=o zU Pn- 

2xf(z) = {az~ l + a- 1 z}f(z)+ 

{-(a + bq- N )z- 1 + a + b + bq' N + a~ l q- N+l - (a^q'^ 1 + b)z}f(qz) + 
bq- N {z- l -q-l + qz}f{q 2 z), 

(31) 

Its solutions in the space spanned by {z k } k v =0 1 are 

2x m = aq m + a- 1 q- m J m = 0, 1, . . . , N - 1, (32) 

f n = v n (ab;q) n (q~ N+1 ;q)n A (?'"; g) fc (a 2 g m ; g) fc (g' m ; g) fc g fc 
JmW 2^ a2n(g . g)n(6a - lg _v + i. g)n Z. (a6;g) fc (g-^;g) fe (g;g) fe 

, , N -^ m (q m - N+1 ; q) k (abq m ; q) k q- mk z k 

[Z,q)m to (ba-iq-^;q) k (q;q) k a^ 

(z; q) m P N ^„ m {za~ 2 q- m ~\ bar x q~ N , a 2 q 2m \q), 

(33) 

where P k (x, a, (3\q) are the little g-Jacobi polynomials. Thus the equation 



(pTD is related to the g-difference equation for these polynomials |14| by the 
transformation P/v_i_ m (x) = f m (xa 2 q m+1 ) / (xa 2 q m+1 ; q) m . 

Finally, consider the dual q-Hahn polynomials. (Other known basic hy- 
pergeometric polynomials leading by the procedure of this section to equa- 
tions of the type (2) can be considered as particular cases of the continuous 
dual q-Hahn or dual q-Hahn polynomials.) These polynomials are defined 
by the recurrence relation (we use N — 1 instead of N in the usual definition 
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0) 



fl>(y)Pn = A n p n+1 + (1 + 7<5g -An - C n )p n + C n Pn-!, 

A n = (1 - q n - N+1 )(l - 7 <f C n = 7 g(l - g»)(<5 - g"""), (34) 
/%) = <T y + 7^<f + \ p-i = 0, n = 0, 1, . . . , TV - 1. 

Setting 7 = abq" 1 , 5 = oft -1 , g _Ar+1 = ac, and multiplying the recurrence 
relation fl34j) by a" 1 , we obtain (|20|) where 2x = t + t" 1 , t = ag y . Thus, 
the first N continuous dual q-Hahn polynomials at ac = q~ N+1 and the dual 
q-Hahn polynomials are the same (up to renaming the parameters). 
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